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Abstract 

We study the autoresonant solution of Dufnng's equation in the 
presence of dissipation. This solution is proved to be an attracting 
set. We evaluate the maximal amplitude of the autoresonant solution 
and the time of transition from autoresonant growth of the amplitude 
to the mode of fast oscillations. Analytical results are illustrated by 
numerical simulations. 

PACS numbers: 02.30.Gp, 02.30.Hq 



Introduction 

By the autoresonance is meant the growth of the amplitude of oscillations 
of a solution to a nonlinear equation under action of an external oscillating 
force. This phenomenon looks like phase locking of a nonlinear oscillator 
through a periodic driver. The phase locking was first suggested to accelerate 
relativistic particles, see (13, [H] . Nowadays the autoresonance is thought of 
as universal phenomenon which occurs in a wide range of oscillating physical 
systems from astronomical to atomic ones [2]. 

In the general case the frequency of a nonlinear oscillator depends on 
the amplitude of oscillations or, what is the same, on their energy. Hence, 
in order to change the energy of a nonlinear oscillator the frequency of the 
external force should be adopted to that of the oscillator. If the force is small 
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then the energy of oscillations changes slowly. In order to remain resonant, 
the frequency of the external force should adopt itself also slowly to the 
frequency of the nonlinear oscillator. Moreover, the backward phenomenon 
proves to occur. Namely, the slow change of the frequency of the driver 
results in that the frequency of the nonlinear oscillator actually follows the 
driver frequency. For a contemporary survey of the mathematical aspects of 
autoresonance we refer the reader to [7j. 

The autoresonance phenomenon in systems with dissipation was earlier 
studied both by means of mathematical models and in physical experiments. 
In particular, the existence of autoresonant solution for the system of three 
coupled oscillators with small dissipation was established in [20J and for the 
system with parametric autoresonance in [5J. In the papers [31 [T2] the thresh- 
old of capture into autoresonance was discussed in the presence of dissipation. 
The resonant phase locking phenomenon in van der Pol Duffmg's equation 
with external driver of slowly varying frequency was studied in [13] . 

In this paper we treat two problems for autoresonance in dissipative sys- 
tems for the time being open. Firstly we prove the existence of an attracting 
set for solution trajectories captured into autoresonance. Such attracting set 
was observed numerically in a number of papers [201 M HU US]- The attrac- 
tor in these systems is a slowly varying steady state solution. The solutions 
captured into autoresonance oscillate around such a solution and lose the en- 
ergy of oscillations because of dissipation. Therefore, all captured solutions 
tend to the steady state solution. Mathematically this means that the slowly 
varying steady state solution is Lyapunov stable. 

The second problem we deal with consists in evaluating the bound of the 
autoresonant growth of solution in the presence of small dissipation in the 
system. Earlier one observed that the amplitude growth of nonlinear oscilla- 
tions in systems with dissipation is bounded, see [TjJl [TU [T5] . /,From physical 
viewpoint the boundedness of autoresonant growth can be easily explained. 
Namely, the work of driver is proportional to the length of trajectory in the 
phase space. If the disspation depends linearly on the velocity then its work 
is proportional to the area described by the phase trajectory. Under the 
growth of energy the area of the phase curve grows faster than its length. It 
follows that even if the dissipation is small, its work exceeds the work of ex- 
ternal force at some moment and the autoresonant growth of solution stops. 
Mathematically this looks like the impossibility of extension of the solution 
under phase capture. What happens is the hard loss of stability and passage 
to fast oscillations. 



2 



This work is aimed at finding an asymptotic expansion for the slowly 
varying steady state solution to the primary resonance equation and at show- 
ing that it is an attracting set for those solutions which are captured into 
the resonance. Moreover, we derive asymptotics for the maximal amplitude 
of oscillations under autoresonance with small dissipation and calculate the 
period of the autoresonant mode in the solution. 

The paper contains seven sections. In the next Section [T] we describe 
mathematical setting of the problem. Section [2] provides a detailed exposition 
of the main results. Section [3] deals with asymptotics for the autoresonant 
mode. In Section H] we discuss the stability of autoresonant growth. In 
Section we study the solution behaviour in the vicinity of the break of 
autoresonant growth. In Section [H] we will look more closely at the break of 
autoresonant growth. Section [7] is concerned with passage from monotone 
autoresonant growth of the amplitude of nonlinear oscillations to fast motion 
solutions. 

1 Setting of the problem 

We study solution of the primary resonance equation 

+ (t-\$\ 2 )$ + i5V = f, (1.1) 

where r is an independent variable, 5 a dissipation parameter and / is a 
parameter related to the amplitude of external force. 

The primary resonance equation is of universal character in mathematical 
description of autoresonance. In the case 5 = it was first introduced in the 
paper [16]. 

The primary resonance equation describes long-term evolution of nonlin- 
ear oscillations under action of a small external force. As but one example 
we mention Duffing's equation with dissipation 

u" + u + bu'-cu 3 = eAcos(ujt) } (1.2) 

where c and A are constants, b and e small positive parameters. The fre- 
quency of oscillations of the right-hand side of the equation depends linearly 
on the time. More precisely, uj = 1 — at which is usually referred to as a 
chirped frequency. The parameter a (called a chirp rate) determines the rate 
of change of the frequency. 
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Duffing's equation (11.21) proves to be the simplest and so the most general 
equation where one observes the phenomenon of autoresonance break because 
of small dissipation. 

For studying autoresonance it is convenient to use the method of two 
scales. The oscillations of the nonlinear equation are observed in the time 
scale t. The amplitude of these oscillations depends in turn on the slow time 
r := e 2 /H. 

The introduction of two time scales enables one to split the evolution of 
solution into two parts, fast and slow ones, using the asymptotic substitution 

u ~ e 1 / 3 !^(r)e i< '*~ r ^ + complex conjugate term 

in equation (11.21) . The standard averaging procedure over the fast time t in 
the leading-order term in e leads to equation (11.11) for the unknown function 
where 5 = e^^b/A, f = A/ (Ay/2) and c = -2y/2. 

This primary resonance equation is often written as the system of equa- 
tions for the amplitude R(t) = \ &(t)\ and the phase <£>(t) = ar g( More 
precisely, 

R! = —5R — / sin <p, 

4 = (r-R 2 )-t cos ^ ( L3 ) 

cf. m. 

The autoresonance or phase locking for the solution of system (I1.3P means 
that ip' = o(l) for r — > oo. This condition along with the second equation 
of system (11.3D determine a low for the amplitude growth which reads R = 
yfr + o(l). The first equation of ( 11. 3ft gives a sufficient condition for the 
instant at which the phase locking is destroyed, namely T* = f 2 /5 2 . Analysis 
of the equation and phase locking condition actually yield an estimate of 
autoresonance growth in a dissipative system with small dissipation, 5 <^ 1, 

see nnnnns]. 

In this paper we construct asymptotics for the slowly varying steady 
state solution of equation (11. 3p with 5 1. We prove that this solution is an 
attracting set for the captured solutions. Moreover, we show asymptotics of 
the maximal value of R and evaluate the instant at which the phase locking 
is destroyed. 

In order to better motivate the problem, we demonstrate results of nu- 
merical simulations for equation (11. ip with 5 > 0. In Figure [1] one can 
observe three stages of evolution for the solution of (II. ip . At the first stage 
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the oscillations are close to some smooth curve. Then, at the second stage 
the solution varies slowly. Finally, at the third stage the solution loses its 
stability and the amplitude of fast oscillations tends to zero. 

The autoresonant growth and rapid decay of the amplitude of & are 
shown in Figure [2j 
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Fig. 1: The solution of (II. ip with parameters / 
condition & \ T =o= 0. 



1, 5 = 0.1 and the initial 



2 Results of the paper 

To formulate the results it is convenient to change both the independent and 
dependent variables by 



e = t5 2 , 



(2.1) 



The equation for ip takes the form 

i*y + (0- \^\ 2 )^ + tS 3 ^ = S 3 f. (2.2) 
Denote R = \ip\ and ip = arg(^) for 9 > and < 8 <C 1. 
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The existence time for the autoresonant mode in the solution of (11. ip is 
evaluated by 

where zq is the first real pole of the Painleve-1 transcendental with zero mon- 
odromy data yi(z, 0,0). Furthermore, the maximal amplitude is estimated 
by 

We are now able to give an explicit description of asymptotics for the 
autoresonant solution which is an attracting set for the solutions captured 
into autoresonance. 

If (f 2 - 9)5~ 1 > 1 then R and y? behave like 



R(0,S) ~ + ^— + 5 5 f2 ~ W 



2e 29^J^9 m 2 (p-ey/ 2 
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86 5 / 2 {f 

ip(9,8) ~ - arctan ( — - ) + 8 — J + 5 2 



J~p~^e~ } 2^6^f] 2 ~^e 8Ve(f 2 - ey/ 2 



jp + 2e)^J 2 ~^~e-2AVe{e-f 2 f 
48(6 -pye 3 / 2 
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To write the asymptotics in a neighbourhood of 9 = f 2 , we change the 
variables by 

V = (9-f 2 )5-\ 
r{r},6) = (R-f)5-\ 
a( V ,6) = ( V -(?>/2)k)5- 1 / 2 . 

The functions r(r], 5) and a(r), S) have the form 



V^V^l J2 V 2 + 4a-l) 
o ■ 



the representations being valid if 8(—l — r/) _1 <C 1. 

Close to 1] = —1 it is convenient to represent the asymptotics in the form 

77 + I 

a = 5 1/2 u(t,5), 



1 Af 2 r-l ... 



In this domain the autoresonant mode of the solution loses its stability The 
leading-order term of u relative to 5 admits the representation 

u{t) ~ 3 3 ' s yi (z,0,0), 

P 1 
r = — z 
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where 0,0) is the Painleve-1 transcendental, see [8], i.e., a special solu- 
tion of the Painleve-1 equation = 6yf + z with asymptotics 



y 1 (z) = ] j^ + 0(z^ 2 ) 
The asymptotic formula for v looks like 

«(r) ~ 
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as 5 — > 0. 

The Painleve-1 transcendental yi(z, 0,0) has poles on the real axis. The 
approximate solution of (11.11) by means of yi(z, 0,0) is valid up to a small 
neighbourhood of the first z of these poles or, what is the same, up to 
r = t := (f 2 /6)z — l/(4/ 2 ). Near the pole the validity domain is determined 
by the inequality 

Vs 

— — < i. 

T - T 

The asymptotics of the solution of (11.11) in a neighbourhood of the pole 
represents by fast non-autoresonant oscillations in the new scale of variable 
6 = f 2 — 5 + 5 2 r + It is convenient to write the unknown functions in 

the form 

R ~ f-± + 5V5p(0, 
3 

if ~ -7T + S(£). 

The function s(£) is a special solution of the equation 

s' = -yjE + f 2 (s -sins), 

such that s — > 0— as ^ — > — oo. The function p(^) is determined from the 
equation 

s' 

Note that the function s(^) depends on a parameter = + (s — sins) 
which tends to as £ — > — oo. 



3 Asymptotics of autoresonant growth 

In this section we construct an asymptotic solution to (12. 2p in the domain 
(f 2 — 0)5~ 1 ^> 1 and 6 > 0. To this end we introduce new unknown functions 
p(0,S) and a(9, 6) related to the amplitude R and phase of the unknown 
function ip by 

R(6,5) = V9 + 5 3 p(9,5), 
<p(9,5) = a(6,5). 
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On substituting these formulas into (12 .2p and separating the real and imagi- 
nary parts of the equation we get 

S 4 p' + + / sin a + 6^-= + 5 3 p = 0, 

2Ve (3.i) 

(6y/0 + 6 4 p)a' + 26p-f cos a + 35 3 Vdp 2 + 5 6 p 3 = 0. 

Assuming 5 to be small, we look for a solution p, a in the form of asymp- 
totic series 

oo 

P (e,8) ~ 

fc ^° (3.2) 

We first derive equations to determine the coefficients of these asymp- 
totic series. For this purpose we substitute (13.21) into equations (13. ip . The 
trigonometric functions in these equations are expanded as Taylor series in 
a neighbourhood of some point cto- Then we equate the coefficients of the 
same powers of parameter 5. As a result we get a recurrent sequence of tri- 
angle systems of linear equations for the unknown coefficients of (13.21) . In 
particular, for the leading-order terms of series we obtain 

\f§ + / sin ao = 0, 
2#Po _ / cos «o = 0, 

which gives p and a . 

On equating the coefficients of 8 we arrive at the system 

cos«o + 1 = 0, 
sin ao = 

which readily yields a\ and p\ by 

1 

«i = 



pi = — 7=(ai 



2V0\ 2V9^/f^6< 
On equating the coefficients of 5 2 we get the system 

—2a 2 cos a + sin a al = 0, 
A9p2 + 2v / #o/ 1 + (cos«o a\ + 2 sin«o 0.2) f = 
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implying 



a 2 = 



P2 = 



On equating the coefficients of 5 3 one still obtains a transparent system for 
two unknown functions a 3 , p 3 

/cosa;o«3 = — Po + ^ ^ cosao a i + 6 sin a a^cci j , 

2V9p 3 - a 3 = —a ' 2 - pi - — cos a p + -a\ = cos ao a^aii 

9 o \/Q 

whose solution is 

1. 3 om 
«3 = -K - 36/) - = ^ ia 2 , 

1 r 'a. ^ V 7 /^ 1 3 (^ 2 + 2)(/ 2 -g) 

and so on. 

Careful analysis of formulas for a k and p k obtained in this way actually 
shows that 

as 9 — > f 2 — 0. From these equalities it follows that the constructed asymp- 
totic expansion is valid for 5(f 2 — 9)^ <C I. 



4 Stability of autoresonant growth 

We will look for a solution which is a partial sum of the asymptotic series 
constructed above, up to remainders p(£,0,6) and a(£,0,5). Namely, we 
consider 

^(0, 5) = (y/9 + Pl {9)5* + p5 3 ) exp (t(a (9) + c*i(0)<S + a5 2 )) (4.1) 
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where < 5 <C 1 and £ = 5~ 3 9 is the fast variable. 

We substitute (14.11) into (12.21) . The task is now to write down the linear 
part of the system for p and a. An easy computation yields the system of 
equations 

Pi = (^P^5-^=)a + m + 0(6% (42) 
a\ = -2V6p + 5a + f 2 (6) + 0(5 2 ). 

The right-hand side of the system has coefficients slowly varying in the 
fast variable £. Solutions of such systems are usually constructed by the 
WKB method, see for instance [IS] . The eigenvalues of the matrix on the 
right-hand side of (14. 2[) are 

A M = ±i V2V(P-9)9 5^ - 5 - + 0(5 3 / 2 ), 

and so the real part of eigenvalues is negative. Hence it follows that the 
asymptotic solution constructed above is stable in linear approximation. Fig- 
ure [3] illustrates this result. 



Irn 5 




Re £ 



Fig. 3: The graph displays the exponential decay of oscillations with zero 
initial data in a neighbourhood of a slowly varying asymptotic solution (bold 
curve) with / = 1 and 5 = 0.05. 
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5 Vicinity of the break of autoresonant 
growth 

We change the variables by 9 = f 2 — Srj. The new independent variable rj is 
stretched with respect to 6. We will look a solution of the form 

p(9,5) = f + 5r( V ,5), 
3 

a(6,8) = -7T + V5a(r],8). 

This substitution leads to the system for two unknown functions r(i], 5) and 
a(rj,5) 

r' = -r + f(cos(V5a)-l)/5, 

a' = 6-^( V -fr-S^)-S^ fs ^ a) . 

v ' f + 5r 

This system can be rewritten in a slightly different form 

fa 2 = 2(-/-r-/ C ° S( ^ Q) ~ 1 + ^ 2/2 , 

5 

2fr = rj-5r*-5^(a'-5-^ fS ™ iV ? a) ). 

V / + or / 

To find asymptotics we substitute formal series in powers of \f5 for a(ri, 5) 
andr(?7,<5). Namely 

oo 

r(r],5) = J2 r ^)$ k/2 , 

oo 

a( V ,5) = J2 a ^v)S k/2 - 



k=0 

On substituting these series into the system we expand both left-hand side 
and right-hand side of the equalities as formal series in powers of \f5. Then 
we equate the coefficients of the same powers of y/d in both series. As a result 
we arrive at a recurrent system of equations for determining the coefficients 
of formal series. For k — it reads 

2/Vo = rj, 
fa 2 = -1-7]. 
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For k — 1 we get T\ — and ai = 0. For = 2 the system is 



2/r 2 



/a a 2 



implying 




and so on. 

The formulas for the coefficients r k and a k are cumbersome. However, 
using the recurrence relations one can see that the coefficients have a singu- 
larity at rj = —1. The greater k, the higher the singularity. This is caused 
by differentiating the square root \J— 1 — 77 and by increasing the nonlinear 
dependence on lower order terms of asymptotics at each step of iteration. 
More precisely, we get 



as 77 — > — 1, provided k — 4 e N. Hence it follows that the constructed series 
is asymptotic for 5 (—1 — i])^ 1 <^ 1. 

6 Break of autoresonant growth 

In a neighbourhood of the point 77 = —1 we change the variables by the 
formula rj = — 1 + tS. The new independent variable r is fast with respect 
to the original variable r]. The solution of the primary resonance equation is 
written in the form 



a k = 0((-l- „)(!-*>/*), 
r k = 0((-l-77)( 4 - fc )/ 2 ) 



a 



T = -777 + 




5 + S 2 v(t,S). 
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Substituting these formulas into the system of equations for a and r, we 
immediately obtain 



v '-U fcos{6u) - f ) + ^-^ +5v = °' 



S2 y V , ,J 2/ g/ 

/ 1 / , r -i -x f sm(5u 

u ~ V v + ^7l( 4 / r - !) + 



We will look for a formal solution to this system in the form of power 
series in 5 

oo 

u(t,5) = J2 u ^)$ k , 

h Z° (6.1) 
v(t,5) = J^^ 6 "- 

k=0 

Substituting these series into the equations and expanding the left-hand sides 
as power series in S, we equate the coefficients of the same powers of 5. This 
leads to a recurrent sequence of differential equations for u k and v%. In 
particular, for u and v we get 

^o + 2M-^(r-^) = 0, 



^0 + ^0+2^-7^) = °- 



2f 2 v 4 / 2 
For ui and v\ the system looks like 

/ 1 1 / 1 \ 2 

«i + 2M = 7^0-^0-^^- — 



u x + -«o«i = ~Vq. 



For m 2 and v 2 the system is 



11 1 

u' 2 + 2fv 2 = jvi-^(l- 4/ 2 t) vq-U!- ^pU , 

l . f f 4 / 2 

U 2 + -Uo«2 = -Ui + — U - -«!, 
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and so on. 

The system of equations for the leading-order terms reduces to the Painleve - 
1 equation. To see this, let 



Then the differentiation of the equation for u leads, by the second equa- 
tion, to the Painleve -1 equations in the standard form y" = 6y 2 + z. Here 
y = y(z, ci, c 2 ) is the first Painleve transcendental, C\ and c 2 are real param- 
eters of the transcendental which are monodromy data, cf. [8], and z is an 
independent variable. 

The solution of the system of equations for the leading-order term is 
determined through the first Painleve transcendental. The parameters of the 
transcendental are defined by making asymptotic expansions consistent. To 
this end we re-expand asymptotic series ( 16. ip in terms of the variable r and 
equate the coefficients of the same powers of 5. Then we get asymptotics of 
the coefficients for r — > — oo, namely 



r + 



4/ 2 






and 
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The asymptotics of the k th correction is 



U2n = 



V2n = 



o(V-r 2n+1 ) 

0(r n+1 ), 



u 2n+ i = 0{V^r^), 
v 2n+ i = 0(r™+ 2 ) 



as r — > — oo. 

The solution u and v with given asymptotics as r — > — oo can be ex- 
pressed through the first Painleve transcendental. The asymptotics of the 
first Painleve transcendental were investigated in [HI El H]- Here it is conve- 
nient to make use of the connection of the asymptotics and the monodromy 
data 



see [8]. Starting with the formula for Uq, one obtains easily an expression for 
Vq(t) from the first equation of the system for uq and Vq- 

We now turn to construction of solutions Uk and t>&. The corresponding 
homogeneous system is 



On differentiating the first equation and substituting V into the second equa- 
tion we arrive at the linearised Painleve -1 equation w" + 2u$w = 0. The 
general solution of this equation is known to be a linear combination of the 
partial derivatives of the first Painleve transcendental in parameters, i.e. 

w = At d ci y(z, ci, c 2 ) + A 2 d C2 y(z, c u c 2 ) 

where A\ and A 2 are arbitrary constants. The asymptotics of y(z, ci, c 2 ) as 
z — >• — oo implies 




U' + 2fV 
V'+ f -u U 



0. 



0. 



(6.2) 



Set 




d Cl y(z,cx,c 2 ) \ 
d C2 y(z,c h c 2 ) \ 



c 1= 
c 2 =0 



c 1= 
c 2 =0 




see 0. 
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The formulas for corrections and t>& can now be obtained by the method 
of variation of constants 



Here <P(t) stands for the fundamental matrix of linearised system (16.21) with 
Wronskian equal to 1. By a is meant an arbitrary real constant satisfying 
a < t := (f 2 /6)z — l/(4/ 2 ), where z is the least real pole of the Painleve 
transcendental. The constants A n and B n are uniquely determined from the 
matching condition for asymptotic solutions. 

The first Painleve transcendental has second order poles on the real axis, 
see [I]. In a neighbourhood of r the constructed asymptotic expansion no 
longer holds. Indeed, we get 

6 



U 

v 



f 2 (r-r r 
6 



/ 3 (r-r )3 
for r close to r . 

The general solution for the first correction u\ and V\ can be represented 
in the form 

"i = TZ fll _ v, + s^J ^^ +O((r-r )- 1 ), 



(r-r ) 3 / 4 (r-r ; 



u i = + T\I + °(( r - r o) )• 

2(r-r ) 4 J b {r-r y 

Here ai is one of the solution parameters. The second independent parameter 
is contained in the smooth part of asymptotics remainder. The parameters of 
the solution are uniquely determined while one constructs it by the method 
of variation of constants. However, in the expansion in a neighbourhood of 
the pole To the parameter at can be included in the pole translation of the 
leading-order term of order 5, namely t\ = r — 5f 2 ai/3. As a result the 
value ro in the expansion of leading-order terms should be replaced by T\ and 
the expansions for u\ and v\ become 



"i = -JI7Z — + ^ T - Tl 



Vl = IH^nf + 0{{T ny2) 
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Thus, the pole of the leading-order term of asymptotics is defined uniquely 
up to 5 2 . More precisely, the pole asymptotics of the perturbed problem is 
determined by singling out summands of order (r — ti)~ 3 in the asymptotics 
of Uk for k > 1. 

The order of singularity at the point t — T\ increases, for the higher order 
corrections depend on lower order corrections in a nonlinear way. For u 2 and 
v 2 we have 
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u 2 ~ 
v 2 ~ 



5/ 6 (r-Ti) 6 ' 
-54 

5/ 7 (r-^i) 7 ' 



One can show that 

U2n-i = 0((T- n )- 2n ), u 2n = Odr-n)-^- 2 ), 
V2n-i = ©((r-n)- 2 "- 1 ), v 2n = 0((r - n)- 4 "- 3 ) 

as r — > ri. 

/From the behaviour of -u^ and t> fe in a neighbourhood of singular point 
we deduce that the constructed asymptotics is valid in the domain 



r - r i 



7 Fast motion 

In a neighbourhood of the singular point T\ the behaviour of the solution 
changes drastically. The solution begins to fastly vary. The new scale of 
independent variable is now 

One introduces new dependent variables S) and s(£, S) by 

P = f-2- + 5V5p(£,5), 

3 f t 7 ' 1 ) 



18 



The genuine independent variable 6 is related to the new independent variable 
£ by the formula 

= / 2 -<5 + <fTi + 

Substituting the expressions for p, a and 9 into original system (13. ip yields 
a system of equations for S and s(£, 5). This system is cumbersome and 
we need not write it here in explicit form. Using the standard procedure of 
perturbation theory we look for the leading-order term of asymptotics in 5 
of the form 

p{£, 8) ~ MO, 

s(£,6) ~ s (0- 
For po(0 an d s o(0 we obtain the system 

Pd + fi 1 ~ cos s ) = 0, 
s' Q + 2fp = 0. 

This system admits the conservation low 

E =p 2 + (sin s - s ). 

Making (17. ip and asymptotics in a neighbourhood of the pole consistent gives 
a condition for p and So, namely both E and So vanish as r] — > — oo. The 
solution of the system for po and so varies fastly and so increases infinitely. 
The variable so stands actually for the argument of the solution in complex 
form 

(/ - 2 j + 5V5p(Z, 6)) exp (3tt + s(£, 5)) . 

Thus, in this mode the phase locking condition fails to hold and the solution 
is not autoresonant. 

Conclusion 

In the paper it is shown for Duffing's equation that in the autoresonance 
domain there exists an asymptotic solution which is stable in linear approxi- 
mation. We evaluate the break time and the maximal amplitude of autores- 
onant solution for the equation with small dissipation. The break of the 
autoresonant mode is accompanied by hard loss of stability and passage to 
fast oscillations. 
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